Detecting synchronization in spatially extended discrete systems by 

complexity measurements 
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The synchronization of two stochastically coupled one-dimensional cellular automata (CA) is 
analyzed. It is shown that the transition to synchronization is characterized by a dramatic increase 
of the statistical complexity of the patterns generated by the difference automaton. This singular 
behavior is verified to be present in several CA rules displaying complex behavior. 

o ■ 
o : 

Despite all the efforts devoted to understand the meaning of complexity, we still do not have an instrument in 
the laboratories specially designed for quantifying this property. Maybe this is not the final objective of all those 
theoretical attempts carried out in the most diverse fields of knowledge in the last years 0, 0> H EL H II] i but, for 
a moment, let us think in that possibility. 

Similarly to any other device, our hypothetical apparatus will have an input and an output. The input could 
be the time evolution of some variables of the system. The instrument records those signals, analyzes them with 
a proper program and finally screens the result in the form of a complexity measurement. This process is repeated 
for several values of the parameters controlling the dynamics of the system. If our interest is focused in the more 
complex configuration of the system we have now the possibility of tuning such an state by regarding the complexity 
plot obtained at the end of this process. 

As a real applicability of this proposal, let us apply it to an a-la-mode problem. The clusterization or synchronization 
of chaotic coupled elements was put in evidence at the beginning of the nineties |jj.llC|. Since then, a lot of publications 
i— h ■ have been devoted to this subject |Tl|. Let us consider one particular of these systems to illuminate our proposal. 

(1) SYSTEM: We take two coupled elementary one dimensional cellular automata displaying complex spatio- 
temporal dynamics |12( . Recently, it has been show that this system can undergo through a synchronization transition 

t-H ' |13| . The transition to full synchronization occurs at a critical value p c of a synchronization parameter p. Briefly the 
numerical experiment is as follows. Two L-cell cellular automata (CA) with the same evolution rule 4> are started 
from different random initial conditions for each automaton. Then, at each time step, the dynamics of the coupled CA 
is governed by the successive application of two evolution operators; the independent evolution of each CA according 
to its corresponding rule <& and the application of a stochastic operator that compares the states a\ and of of all the 
\ cells, i = 1, ...L, in each automaton. If a\ = a^, both states are kept invariant. If a\ ^ of, they are left unchanged 
. with probability 1 — p, but both states are updated either to a} or to of with equal probability p/2. It is shown in 
reference [l3( that there exists a critical value of the synchronization parameter (p c = 0.193 for the rule 18) above 
which full synchronization is achieved. 

(2) DEVICE: We choose a particular instrument to perform our measurements, that capable of displaying the value 
of the LMC complexity (C) |LJ|. The statistical complexity C is defined as follows, 
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where {pi} represents the set of probabilities of the N accessible discrete states of the system, with pi > , i = 1, ■ ■ ■ , iV, 
and k is a constant. If k = 1/logN then we have the normalized complexity. C is a statistical measure of complexity 
that identifies the entropy or information stored in a system and its disequilibrium, i.e., the distance from its actual 
state to the probability distribution of equilibrium, as the two basic ingredients for calculating the complexity of 
a system. This quantity vanishes both for completely ordered and for completely random systems giving then the 
correct asymptotic properties required for a such well-behaved measure. The calculation of C has been useful to 
successfully discern many situations regarded as complex in discrete systems out of equilibrium 0, 0, 0, fisl IT^ | . 

(3) INPUT: In particular the evolution of two coupled CA evolving under the rules 22, 30, 90 and 110 is analyzed. 
The pattern of the difference automaton will be the input of our device. In Fig. 1 it is shown it for a coupling 
probability p = 0.23, just above the synchronization transition. The left and the right plots show 250 successive states 
of the two automata, whereas the central plot displays the corresponding difference automaton. Such automaton is 
constructed by comparing one by one all the sites (L = 100) of both automata and putting zero when the states a} 
and of, i = 1 . . . L, are equal or putting one otherwise. It is worth to observe that the difference automaton shows 
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an interesting complex structure close to the synchronization transition. This complex pattern is only found in this 
region of parameter space. When the system if fully synchronized the difference automaton is composed by zeros in 
all the sites, while when there is no synchronization at all the structure of the difference automaton is completely 
random. 

(4) METHOD OF MEASUREMENT: How to perform the measurement of C for such two-dimensional patterns 
has been presented recently in Ref. [23] ■ We let the system evolve until the asymptotic regime is attained. The 
variable af in each cell of the difference pattern is successively translated to an unique binary sequence when the 
variable i covers the spatial dimension of the lattice, i = 1, . . . , L, and the time variable n is consecutively increased. 
This binary string is analyzed in blocks of n Q bits, where n can be considered the scale of observation. The accessible 
states to the system among the 2™° possible states is found as well as their probabilities. Then, the magnitudes H, 
D and C are directly calculated and screened by the device. 

(5) OUTPUT: The results of the measurement are shown in Fig. 2. The normalized complexity C as a function 
of the synchronization parameter p is plotted for different coupled one-dimensional CA that evolve under the rules 
22, 30, 90 and 110 , which are known to generate complex patterns. All the plots of Fig. 2 were obtained using 
the following parameters: number of cell of the automata, L — 1000; total evolution time, T = 600 steps. For all 
the cases and scales analyzed, the statistical complexity C shows a dramatic increase close to the synchronization 
transition. It reflects the complex structure of the difference automaton and the capability of the measurement device 
here proposed for clearly signaling the synchronization transition of two coupled CA. 

These results are in agreement with the measurements of C performed in the patterns generated by a one-dimensional 
logistic coupled map lattice in Ref. 20] . There the LMC statistical complexity (C) also shows a singular behavior 
close to the two edges of an absorbent region where the lattice displays spatio-temporal intermittency. Hence, in 
our present case, the synchronization region of the coupled systems can be interpreted as an absorbent region of the 
difference system. In fact, the highest complexity is reached on the border of this region for p w 0.2. The parallelism 
between both systems is therefore complete. 

Finally, let us remark that the critical parameter for switching synchronization in these systems is similar for all of 
them (p c w 0.2). This fact could deserve some operational explanation in a future work. 
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FIG. 1: Spatio-temporal patterns just above the synchronization transition. The left and the right plots show 250 successive 
states of the two coupled automata and the central plot is the corresponding difference automaton for the rules 22, 30, 90 and 
110. The number of sites is L = 100 and the coupling probability is p = 0.23. 
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FIG. 2: The normalized complexity C versus the coupling probability p for different scales of observation: n = 1 (°),n = 4 
(*) and n = 6 (□). C has been calculated over the last 300 iterations of a running of 600 of them for a lattice with L — 1000 
sites. The synchronization transition is clearly depicted around p » 0.2 for the different rules. 



